Testing replica predictions in experiments 
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Abstract 

We review the predictions of the replica approach both for the statics and for the 
off-equilibrium dynamics. We stress the importance of the Cugliandolo-Kurchan off- 
equilibrium fluctuation-dissipation relation in providing a bridge between the statics 
and the dynamics. We present numerical evidence for the correctness of these relations. 
This approach allows an experimental determination of the basic parameters of the 
replica theory. 

When we suddenly decrease the temperature in an Hamiltonian system, many interesting 
phenomena happen if the initial and the final temperatures correspond to different phases. 
When the low temperature phase can be characterized by a simple order parameter (e.g. the 
magnetization for ferromagnets) we find the familiar phenomenon of spinodal decomposition 
characterized by growing clusters of different phases. There is a dynamical correlation length 
(i.e. the size of the clusters), which increases as a power of the time t after the quench, and 
the energy approaches equilibrium with power like corrections (e.g. E(t) ~ + At~ l l 2 ). 
In this region aging phenomena are also present [|TJ . 

The situation is more intriguing in the case of structural glasses and spin glasses where 
the low temperature phase cannot be characterized in term of a simple order parameter. 
Remarkable progresses in understanding the off-equilibrium dynamics and its relations to 
the equilibrium properties has been done by noticing that a crucial off-equilibrium feature is 
the presence of deviations from the well known equilibrium fluctuation-dissipation relations. 
On the basis of analytic results for soluble models it has been conjectured that we can define 
a function X(C), C being an autocorrelation function at different times 0, [|, [|. This 
function characterizes the violations of the fluctuation-dissipation theorem (which is correct 
only at equilibrium). It is remarkable that (at least in the case of spin glasses) the function 
X(C) is equal to the function x(q) (q being the overlap of two spin configurations) which 
plays a central role in the equilibrium computation of the free energy ||. 

Let us first recall the description of these glassy systems at equilibrium according to 
the predictions of the replica theory. We consider a systems and we denote by C a generic 
configuration of the system. For simplicity we will assume that there are no symmetry in 
the Hamiltonian, in presence of symmetries the arguments must be slightly modified. It is 
useful to introduce an overlap q(C,C). There are many ways in which an overlap can be 
defined; for example in spin system we could define 

q= N > (1) 



1 



N being the total number of spins or particles and a and r are the two spin configurations. 
In a liquid a possibility is given by 



Ei=i,jv Efc=i,jv f{x(i) - y{k)) 
q= N ' (2) 

where / is a function which decays in a fast way at large distances and is substantially 
different from zero only at distances smaller that the interatomic distance (x and y are the 
two configurations of the system). 

In the high temperature phase for very large values of N the probability distribution of 
the overlap (P^-(q)) is given by 

P N (q)*6(q-q*). (3) 

In the low temperature phase -P/v(<z) depends on N (and on the quenched disorder, if it is 
present). When we average over N we find a function P(q) which is not a simple delta 
function. In all known case one finds that 

P{q) = a m S(q - q m ) + a M 5(q - q M ) +p{q), (4) 

where the function p(q) does not contains delta function and its support is in the interval 
[q m ,qM]- 

The non triviality of the function P(q) (i.e. the fact that P(q) is not a single delta 
function and consequently q is an intensive fluctuating quantity) is related to the existence 
of many different equilibrium states. Moreover the function P/v(g) changes with N and its 
statistical properties (i.e the probability of getting a given function P N (q)) can be analytically 
computed |5|, |J . 

In this equilibrium description a crucial role is given by the function x(q) defined as 

x(q) = f P{q')dq' . (5) 

In the simplest case the function p(q) is equal to zero, i.e the function P(q) has only two 
delta functions without the smooth part. In this case, which correspond to one step replica 
symmetry breaking, there are many equilibrium states, labeled by a, and the overlaps among 
two generic configurations of the same state and of two different states are respectively q^ 
and q m . The probability of finding a state with total free energy / is proportional to 

exp(m/3(/ -/*)), (6) 

where /r is a reference free energy and m is the value of x(q) in the interval [q m , qi^]- 

In the more complicated situation where the function p(q) is non zero, couples of different 
states may have different values of the overlaps. The conjoint probability distribution of the 
states and of the overlaps can be described by formulae similar to eq. @, but more complex 
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Although these predictions are quite clear, it is not so simple to test them for many 
reasons: 

• They are valid at thermal equilibrium, a condition that is very difficult to reach for 
this kind of systems. 

• Experimentally is extremely difficult to measure the values of the microscopic variables, 
i.e all the spins of the system at a given moment. These measurements can be done 
only in numerical simulations, where the observation time cannot be very large. 



A very important progress has been done when it was discovered || that the function X, 
which describes the violations of the fluctuation dissipation theorem, is equal to the function 
x which is relevant for the statics. This equality is very interesting because function X(C) 
can be measured relatively easily in off-equilibrium simulations 0. 

The temperature dependence of the function X(C) (or equivalently x(q)) is interesting 
also because rather different systems can be classified in the same universality class according 
to the behaviour of this function. It has been conjectured long time ago that the equilibrium 
properties of glasses are in the same universality class of some simple generalized spin glass 
models §, |. 

Let us be more precise. We concentrate our attention on a quantity A(t). We suppose 
that the system starts at time t = from an initial condition and subsequently it remains at 
a fixed temperature T. If the initial configuration is at equilibrium at a temperature T' > T, 
we observe an off-equilibrium behaviour. We can define a correlation function 



C(t,t w ) = {A(t w )A(t + t w )) (7) 

and the response function 

filA(+J-i \\ 

(8) 



e=0 



5e(t w ) 

where we are considering the evolution in presence of a time dependent Hamiltonian in which 
we have added the term / dte(t)A(t). 

The usual equilibrium fluctuation-dissipation theorem (FDT) tells us that 

G-M = -f££>, (9, 

where 

G eq (t) = lim G(t,t w ), C eq (t) = lim C(t,t w ). (10) 
It is convenient to define the integrated response: 

R{t,t w )= C drG(t-T,t w + T), R eq {t) = lim R{t,t w ), (11) 

JO t m ^co 

R(t, t w ) is the response of the system at time t + t w to a field acting for a time t starting at 
t w . The usual FDT relation becomes 

R e "(t) = (3{C eq {t) - C eq (0)). (12) 

The off-equilibrium fluctuation-dissipation relation states that the response function 
and the correlation function satisfy the following relation for large t w : 

rC(0,t w ) 

R(t,t w ) w (3 / X(C)dC. (13) 

Jc(t,t w ) 

If we plot R versus (3C for large t w the data collapse on the same universal curve and the 
slope of that curve is —X(C). The function X{C) is system dependent and its form tells us 
interesting information, in the case of three dimensional spin glasses. 
We must distinguish two regions: 

• A short time region where X(C) = 1 (the so called FDT region) and C belongs to the 
interval J (i.e. C\ < C < C 2 X 
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Figure 1: The response R times T versus C at T = 0.7 for the three dimensional Ising spin 
glass [10]. The curve is the prediction for function R(C) obtained from the equilibrium data. 
The straight line is the FDT prediction. We have plotted the data of the two runs: t w = 10 5 , 
and £.„ = 10 4 . 



A large time region (usually t = 0(t w )) where C I and X{C) < 1. In the same 
region the correlation function often satisfies an aging relation i.e C(t, t w ) depends only 
on the ration s = t/t w in the region where both t and t w are large: C(t, t w ) ~ C a {t/t w ). 



In the simplest non trivial case, i.e. one step replica symmetry breaking || |9j 
function X{C) is piecewise constant, i.e. 



X(C) 



m 



for C El, X{C) = 1 for C £ I. 



the 



(14) 



One step replica symmetry breaking for glasses has been conjectured in ref. [[|. 

In all known cases in which one step replica symmetry holds, the quantity m vanishes 
linearly with the temperature at small temperatures. It often happens that m = 1 at T = T c 
and m(T) is roughly linear in the whole temperature range. 

Let us consider the case of spin glasses at zero magnetic field (in this case the replica 
symmetry is fully broken |I0|). The natural variable to consider is a single spin (A = <jj)). 
In this case the correlation C(t, t w ) is equal to the overlap among two configurations at time 
t and t w : 

J2i=N CT i('t) 'i{'tw) 



C(t, t v 



N 



(15) 



The response function is just the magnetization in presence of an infinitesimal magnetic field. 
In this case the situation is quite good because there are reliable simulations for the system 
at equilibrium jL0|. 

In fig. (1) (taken from ]TT[) we plot the prediction for the function R versus C, obtained 
at equilibrium (i.e. using the equilibrium probability distribution of the overlaps, P(q)) by 
means of a simulation of a 16 3 lattice using parallel tempering |T2|, [TD[. The simulation has 
been done with the help of the APE100 supercomputer |U| and involves the study of 900 
samples of a L = 16 lattice. 
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During the off-equilibrium simulations ]TT[ m a first run without magnetic field the au- 



tocorrelation function has been computed. In a second second run from t — until t = t w 
the magnetic field is zero and then (for t > t w ) there is an uniform magnetic field of small 
strength ho- The starting configurations were always chosen at random (i.e. the system is 
suddenly quenched from T = oo to the simulation temperature T). 

In fig. (1) there are the results of the off-equilibrium simulations Jll[ where t w = 10 5 



and t w = 10 4 , with a maximum time of 5 • 10 6 Monte Carlo sweeps. The lattice size in was 
64, and T = 0.7 (well inside the spin glass phase, the critical temperature is close to 1.0). 
We plot the response function R times T (in this case R is equal to m/ho) against C(t,t w ). 
We have plotted also a straight line with slope —1 in order to control where the FDT is 
satisfied. Finally we have plotted two points, in the left of the figure, that are obtained with 
the infinite time extrapolation of the magnetization. 

The agreement among the absolute theoretical predictions (no free parameters) coming 
from the statics and the dynamical numerical data is quite remarkable. These data show 
the correctness of the identification of the functions x of the statics and X of the dynamics. 

Let us now go to the case of glass forming materials. I will present the data for binary 
mixture of soft spheres |]14| . Theoretically there have been many speculations the glass tran- 



sition is described by one step replica symmetry breaking [§, ||, [16|, |TJ]. Here the equilibrium 
properties are no so well known as in spin glasses, although there are some evidence that 
the homologous of the function P(q) is non trivial |J. On the other side, as we shall see, 
off-equilibrium simulations |T5| show that the function X(C) seems to be given by the one 
step formula fll~4]) with an approximate linear dependence of m on the temperature. 

We consider a mixture of soft particles of different sizes. Half of the particles are of type 
A, half of type B and the interaction among the particles is given by the Hamiltonian: 

/Mo+awy^ (16) 



i<k 



where the radius (a) depends on the type of particles. This model has been carefully studied 
in the past jL4], |j~5| . The choice ub/^a — 1-2 strongly inhibits crystallisation and the system 
goes into a glassy phase when it is cooled. Using the same conventions of the previous 
investigators we consider particles of average radius 1 at unit density. It is usual to introduce 
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Figure 3: The quantity m = as t w = 2048 as function of the temperature. The straight 
line is the prediction of the approximation m(T) = T/T c . 



the quantity T = (3 A . For quenching from T = oo the glass transition is known to happen 
around T r = 1.45 |14 . 



The best quantity we can measure to evidenziate off-equilibrium effects is the diffusion 
of the particles: 

A(t, u s S^M^fa + '>!'>, (17) 

The usual diffusion constant is given by D = lim^oo A(t, t w )/t. 

The other quantity we measure is the response to a force. At time t w we add to the 
Hamiltonian the term e f • where / is vector of squared length equal to d — 3 and we 
measure the response 



R(t w , t) 



d{f-x k {t w + t)) ( 



de 



(f • x k (t w + t)) t 
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for sufficiently small e. The usual fluctuation theorem tells that at equilibrium /3A eq (t) = 
R eq {t). 

In the following we will look for the validity in the low temperature region of the gener- 
alized relation f3X(A) = dR/dA. This relation (with X ^ 1) can be valid only in the region 
where the diffusion constant D is equal to zero. Strictly speaking also in the glassy region 
D / 0, because diffusion may always happens by interchanging two nearby particles (D is 
different from zero also in a crystal); however if the times are not too large the value of D 
is so small in the glassy phase that this process may be neglected in a first approximation. 

The simulations we present are done using a Monte Carlo algorithm, which is a discretized 
form of a Langevin dynamics. In fig. ^|we show R versus (3 A at t w = 2048 and t w = 8192 for 
T = 1.6 and t < At w at iV = 66. We also show the data for t w = 2048 at iV = 130. We do not 
observe any significant systematic shift in this plot among three data sets. We distinguish 
two linear regions with different slope as expected from one step replica symmetry breaking. 
The slope in the first region is compatible with 1, as expected from the FDT theorem, while 
the slope in the second region is near 0.62. Also the data at different temperatures for all 
values of T > 1.5 show a similar behaviour. The value of R, in the region where the FDT 
relation does not hold, can be very well fitted by a linear function of A as can be seen in fig. 

The region where a linear fit (with m < 1) is quite good corresponds to t/t w > 0.2. 

The fitted value of m = dR/d((3A) is displayed in fig. (3). When m becomes equal to 
1, the fluctuation-dissipation theorem holds in the whole region and this is what happens at 
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higher temperatures. The straight line is the prediction of the approximation m(T) = T/T c , 
using T c = 1.45. 

All the results are in very good agreement with the theoretical expectations based on 
our knowledge extracted from the mean field theory for generalized spin glass models. The 
approximation m(T) = T/T c seems to work with an embarrassing precision. We can conclude 
that the ideas developed for generalized spin glasses have a much wider range of application 
than the models from which they have been extracted. It likely that they reflect quite 
general properties of the phase space and therefore they can be applied in cases which are 
very different from the original ones. 

The most interesting development would be to measure experimentally the function X 
both in spin glasses and in structural glasses. Clearly the most difficult task is the measure- 
ment of the fluctuations. In spin glasses it is clear how it should be done: the measurement 
of the thermal fluctuations of the magnetization is a delicate, but feasible experiment. In 
the case of structural glasses some ingenuity is needed in planning the experiments. A open 
interesting possibility would to do the measurements in the case of rubber, where a transition 
with similar characteristics should take place. 
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